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ABSTRACT 

We present the detection of the reactive ion CO"*" towards the prototypical 
starburst galaxy M 82. This is the first secure detection of this short-lived ion in 
an external galaxy. Values of [CO"'"]/[HCO^]>0.04 are measured across the inner 
650 pc of the nuclear disk of M 82. Such high values of the [CO+]/[HCO"'"] ratio 
had only been previously measured towards the atomic peak in the reflection 
nebula NGC 7023. This detection corroborates that the molecular gas reservoir 
in the M 82 disk is heavily affected by the UV radiation from the recently formed 
stars. Comparing the column densities measured in M 82 with those found in 
prototypical Galactic photon-dominated regions (PDRs), we need ~20 clouds 
along the line of sight to explain our observations. We have completed our model 
of the molecular gas chemistry in the M 82 nucleus. Our PDR chemical model 
successfully explains the [CO^]/[HCO"'"] ratios measured in the M 82 nucleus 
but fails by one order of magnitude to explain the large measured CO^ column 
densities (~l-4xl0^^ cm~^). We explore possible routes to reconcile the chemical 
model and the observations. 

Subject headings: galaxies: individual (M 82) — galaxies: nuclei — galaxies: 
starburst — ISM: molecules — ISM: abundances - radio lines: galaxies 



^Observatorio Astronomico Nacional (OAN), Apdo. 112, E-28803 Alcala de Henares (Madrid), Spain 

^Laboratoire d'Etude du Rayonnement et de la Matiere en la Astrophysique, UMR 8112, CNRS, Ecole 
Normale Superieure et Observatoire de Paris, 24 rue Lhomond, 75231 Paris Cedex 05, France 

■^Departamento de Fisica, UEM, Urb. El Bosque, E-28670 Villaviciosa de Odon (Madrid), Spain 

''^Herschel Science Centre, European Space Astronomy Centre (ESAC), Urb. Villafranca del Castillo P.O. 
Box 50727, E-28080 Madrid, Spain 

^LUTH and UMR 8102 du CNRS, Observatoire de Paris, 92195 Meudon Cedex, France 



- 2 - 



1. Introduction 

M 82 is one of the nearest and brightest starburst galaxies. Located at a distance of 
3.9 Mpc, and with a bolometric luminosity of 3.7x10^° L0, it has been extensively studied 
in many molecules. Several studies reveal that the starburst has heavily influenced the inter- 
stellar medium in M 82 by producing high UV and cosmic ray fluxes. Recent interferometric 
observations at miUimeter wavelengths (Garcia-Burillo et al. 2001) show that while the chem- 
istry of the molecular gas in the disk-halo interface is dominated by shocks, the chemistry of 
the molecular gas in the M 82 disk seems to be dominated by the intense UV flux. Garcfa- 
Burillo et al. (2002) obtained a high-angular-resolution image showing widespread enhanced 
HCO abundance ([HCO]/[H^'^CO^]~3.6) across the whole M 82 disk which was interpreted 
in terms of a giant PDR of 650 pc size. Fuente et al. 2005 (hereafter Paper I) observed a 
selected set of PDR tracers (CN,C2H,H0C~'" and C-C3H2) in three positions across the M 82 
disk. They measured [CN]/[HCN]~5 in the inner 650 pc galaxy disk. This large value of the 
[CN]/[HCN] ratio is only reached in the most heavily exposed layers of a PDR to the UV 
radiation (Fuente et al. 1993; Fuente & Martm-Pintado 1997; Sternberg & Dalgarno 1995; 
Boger & Sternberg 2005). Furthermore, we detected the HOC"*" 1— >0 line and obtained an 
[HCO+]/[HOC+] ratio of ~ 40. Such a low [HCO+]/[HOC+] ratio is only expected in highly 
ionized molecular gas (X(e~)>10~^) either by UV photons (PDRs) and/or X-rays (XDRs) 
(Fuente et al. 2003; Rizzo et al. 2003; Usero et al. 2004). 

In this paper we report the CO"*" detection in the nucleus of M 82. This is the first 
secure C0+ detection in an extragalactic object since previously it has only been tentatively 
detected towards the active galactic nucleus (AGN) Cygnus A (Fuente et al. 2000). CO"'" 
is a reactive ion that can only survive in the highly ionized layers of photon-dominated 
and X-rays dominated regions (PDRs and XDRs) (see e.g. Sternberg & Dalgarno 1995). 
In our Galaxy, CO"*" has only been detected in a handful of objects which are well known 
prototypical PDRs (NGC 7027 and M17 SW: Latter et al. 1993, Fuente et al. 2003; the 
Orion Bar: Storzer et al. 1995, Fuente et al. 2003; NGC 7023: Fuente & Martin-Pintado 
1997, Fuente et al. 2003; Mon R2 and G29.96-0.02: Rizzo et al. 2003; S140 and NGC 2023: 
Savage & Ziurys 2004). In fact, CO^ may be the best molecular tracer of the outermost layers 
(A^,>2 mag) of PDRs. In contrast with other molecular PDR tracers like CN and HOC"*", 
C0+ is exclusively formed in these layers by reactions of C"*" with OH. The formation of the 
chemically related ion HOC"'" is favored in this region but can also be formed at a smaller 
rate in molecular clouds. In addition to the CO"'" detection, we present observations of the 
high excitation HCO+ 3^2, HOC+ 3^2 and the CH3OH 5k,k'^4:k,k' lines. 
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2. Observations and Results 

The observations were carried out in June and November 2004 with the IRAM 30m 
radiotelescope at Pico de Veleta (Spain). We used 2 SIS receivers tuned in single-sideband 
mode in the 1 mm band. The observed transitions are: HCO+ 3— >2, HOC+ 3^2, CO"'" 
N=2^1 J=5/2^3/2 and J=3/2^1/2, and CH3OH 52^42 and 5_2^4_2 E . In Fig. 1 we 
present the observed spectra and the Gaussian fits are shown in Table 1. The intensity scale 
is main beam brightness temperature. The forward efficiency (77//), main-beam efficiency 
(Vmb) and half-power beam width (HPBW) of the telescope are 0.91,0.52 and 11" at 235 
GHz and 0.88, 0.46 and 9" at 260 GHz. Pointing was checked every two hours by observing 
a reference source. We observed three positions across the M 82 disk in the CO"*" lines: the 
nucleus (RA(2000): 09''55"51.9^ Dec(2000): 69°04'47.11") (hereafter referred to as Center) 
and the two peaks in the HCO emission [offsets (-1-14", -|-5") and (—14",— 5") hereafter referred 
to as £'and Irrespectively]. Only E'and Center were observed in the HCO"*" 3— >2 and HOC+ 
3— >-2 hnes. The HOC+ 3— >2 line has been detected towards E. This detection constitutes 
a further corroboration of the HOC"*" 1— >^0 detection reported in Paper I. The CH3OH hne 
was only observed towards Center. 



2.1. CO+ detection 

CO"*" has a ground electronic state in which each rotational level is split into two 
fine-structure levels with J=N±l/2. The N=l— >0 rotational line is heavily obscured by 
the O2 hne at 118 GHz and cannot be observed from ground-based telescopes. The most 
intense transitions of the N=2^1 rotational spectrum are J— 5/2^3/2 at 236.062 GHz and 

J=3/2-^l/2 at 235.789 GHz. Since they are very close in frequency, the two lines can be 
observed simultaneously. In the optically thin limit, the intensity ratio I(236.062)/I(235.789) 
is 1.8. The detection of the two lines with the expected intensity ratio provides supporting 
evidence of the reality of the CO"*" detection. 

The C0+ N=2^1 J=5/2-^3/2 line has been detected towards E (>4o-) and Center 
(>8cr) and very tentatively (~3cr) towards W. Furthermore, we have a ~3o" detection of 
the weakest component towards Center. The 236.062 GHz CO"*" line is blended with the 
5_2,4^4_2,3 and 52,3^42,2 E lines of ^3CH30H (at 236.062 and 236.063GHz respectively). To 
confirm the CO"*" detection and estimate the possible contamination of the 13-methanol lines, 
we have observed towards Center the same transitions of the abundant isotope ^^CHsOH 
and obtained Tmb~10 mK. Assuming a ^^C/^^C ratio >50 (Mao et al. 2000), the intensity 
of the CH3OH line should be <0.2 mK towards Center, i.e., at least a factor of 30 lower than 
the intensity observed in the C0+ spectrum. Therefore, we can conclude that the detected 
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emission at 236.062 GHz corresponds to the J=3/2— >l/2 line of 00+. 

2.2. Column density ratios 

The physical conditions of the molecular gas are estimated by fitting the intensities of the 
H^^CO''" 1— >^0 and HCO+ 3— >2 lines using an LVG code. For these calculations, we assume 
Tfc=50 K (WeiB et al. 2001) and [HCO+]/[Hi3CO+]=89. The ^'^C/^^C ratio in the M 82 
nucleus is not well known. Mao et al. (2000) derived a ^^C/^^C ratio between 50 and 75 from 
multiline CO observations using an LVG code. However they obtained a larger value of the 
i2q^i3q YoXiio when applying a PDR model to the same data. For consistency with Paper 
I, wc adopt the canonical value [HCO^]/[H^'^CO^]=89. The derived molecular hydrogen 
densities arc 3-8x10^ cm~^. These densities arc in agreement with those derived from the 
CN and HCN lines in Paper I within the uncertainties inherent to this kind of calculations. 
Assuming these densities and an emission size of 6" for E, W and Center we obtain the 
H^^GO+ and HOG+ column densities shown in Table 2. The size of the emission towards E, 
VFand Center has been derived from the interferometric HCO and H^^CO+ images published 
by Garcia-Burillo et al. (2002) (see Fig. 1). To calculate the GO"*" column density, we assume 
optically thin emission and use the Local Thermodynamic Equilibrium (LTE) approximation 
with Tj.ot=10 K. This rotation temperature has been estimated using a LVG code applied 
to a linear molecule with the same dipole moment as C0+ and n(H2) = 10^ cm~^. Assuming 
a size of 6" for the 00"*" emission we derive CO"'" column densities ~l-4xl0^^ cm~^ across 
the M 82 disk. 

We can calculate the [HCO+]/[HOC"^] ratio towards E using the HCO"*" 3— s>2 and 
HOC"*" 3-^2 lines. Since both lines lie at the same frequency and have similar dipole moments, 
[HCO+]/[HOC+]«I(HCO+ 3^2)/I(HOC+ 3^2) ratio. We obtain [HCO+]/[HOC+]fti48 in 
excellent agreement with our previous estimates (Paper I). 

In Table 2 we compare the column density ratios in M 82 with those measured in some 
prototypical Galactic PDRs. The [CO"*"] /[HCO"*"] ratio is larger than 0.04 all across the 
M 82 disk. This is one of the largest values of the [CO"'"]/[HCO+] ratio measured thus far. 
Values of the [CO+j/fHCO"*"] ratio larger than 0.01 are only found towards the atomic peaks 
in prototypical Galactic PDRs, and such a high value of the [GO''"]/[HGO"'"] ratio has only 
been measured in the most exposed layers of the PDR associated with the reflection nebula 
NGG 7023. This result is consistent with the values of the [HCO+]/[HOC+] and [CN]/[HCN] 
ratios previously measured in the M 82 disk (Paper I). We estimated [HCO"'"]/[HOG"'"]~40 
and [CN]/[HCN]>5 across the M 82 nucleus. These values are also similar to those measured 
in the PDR peak towards NGG 7023. The very favorable geometry of the PDR associated 
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with this reflection nebula allowed us to detect the outermost layers of the PDR (Fuente 
et al. 1993, 1996a). Our observations suggest that the bulk of the dense molecular gas in 
M 82 is surviving in a similar environment to that found in the HI/H2 transition layer of 
this photon-dominated region. 

We can also compare the C0+, HOC+ and CN column densities in M 82 with those 
derived in Galactic PDRs, although the derived column densities arc more uncertain than 
the column density ratios since the former depend on the assumed beam filling factor. In 
Galactic PDRs, the C0+ column densities are quite uniform taking a value of ~10^^ cm~^ 
in all star forming regions regardless of the incident UV field in a range of 3 orders of 
magnitude. Furthermore, the [CO^] /[HOC'^] ratio is ~0.5-9 (Rizzo et al. 2003; Savage & 
Ziurys 2004). The [CO+]/[HOC+] ratio measured in M 82 is similar to those found in the 
PDRs associated with star forming regions in our Galaxy (see Table 2). This suggests a 
similar CO+/HCO+/HOC"'" chemistry and a similar origin for the reactive ions. Since the 
CO"'" column density is a factor 20-40 larger in M 82 than in Galactic PDRs, we need about 
20-40 PDRs along the line of sight to account for our observations. In the scenario of clouds 
immersed in a pervasive UV field, this implies about 10-20 clouds along the line of sight 
which is a reasonable number for an edge-on galaxy. 



3. Chemical Model 

To have a deeper insight into the physics and chemistry of the molecular clouds in 
M 82, we have modeled their chemistry using the PDR Le Bourlot's semi-infinite plane par- 
allel model (Le Bourlot et al. 1993) and the same physical conditions as in Paper I: Go=10^ 
in units of Habing field, n=nj/+2xn//2=4x 10^ cm~^ and a cosmic ray fiux of 4xl0~^^ s~^. 
As discussed in Paper I, this model explains with reasonable success all the molecular abun- 
dances observed in M 82 and should therefore account for our C0+ detection. 

The model results for C0+ and HCO+ arc shown in Fig. 2. The N(CO+)/N(HCO+) 
ratio (hereafter rco+) is very high (>0.1) for A^,<0.5 mag. Then, rco+ remains constant 
and equal to ~0.1 until A^<5 mag. This is because in this range of visual extinction, 
0.5<A„<4.5 mag, rco+ is determined by the C0+ and HCO"*" abundances in the most 
external layers (A^<1 mag) of the PDR. For higher extinctions, rco+ decreases because of 
the rapid increase of the HCO"*" abundance. In our plane-parallel model, the values of rco+ 
observed across the M 82 nucleus (rco+>0.04) are found only for A^<6.5 mag (see Fig. 2). 
In an external galaxy, one does not expect to have a single PDR but a population of clouds 
(or cloudlets) immersed in an intense UV field (Paper I, Garcia-Burillo et al. 2002). In 
this scenario, our model results imply that the individual cloudlets have Ntoj~1.3 10^^ cm~^. 
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Thus, our 00+ observations corroborate the scenario for the M 82 nucleus proposed in 
Paper I of a highly fragmented interstellar medium in which the dense cores (n~4xl0^ cm~^, 
Ntot~1.3xl0^^ cm~^ ) are bathed by an intense UV field (Go=10^ Habing fields). 

Thus far, we have only compared the observed and model predicted molecular column 
density ratios. We can also compare the molecular column densities. The predicted CO"*" 
column density is ~3xl0^'^ cm^^ for A,u=6.5 mag. This value is a factor of 20-40 lower 
than the CO"'" column densities observed in the prototypical Galactic PDRs. Furthermore 
it is 3 orders of magnitude lower than the C0+ column densities measured in M 82 and 
an unrealistic large number of PDRs along the line of sight would be required to explain 
our C0+ observations. To have a deeper insight into the cause of this discrepancy between 
theoretical predictions and the observations, we have also compared the predicted CN column 
densities with the observed ones (see Table 1). Like C0+, CN is a good tracer of PDRs and 
can be used to estimate the number of PDRs along the line of sight. Our model predicts 
N(CN)~1.5xlO''"^ cm~^ for Ai,~6.5 mag. This value agrees within a factor of 2 with those 
observed in prototypical Galactic PDRs like the Orion Bar and NGC 7023. Comparing with 
the CN column densities observed in the M 82 nucleus, we need about 20-40 individual 
cloudlets along the line of sight to account for our observations (see also Boger & Sternberg 
2005). This is a reasonable number of cloudlets for an edge-on galaxy. Thus, there is a 
reasonable agreement between model predictions and observations for CN in both. Galactic 
and extra-galactic PDRs. However, the chemical model falls short by more than one order 
of magnitude of accounting for the CO"*" column densities measured in Galactic PDRs and 
the M 82 nucleus. 

The failure of chemical models to account for the observed reactive ions column densities 
is a long standing problem (Black 1998; Fuente et al. 2000). The chemistry of reactive ions 
is very sensitive to the gas physical conditions in the HI/H2 transition layer. In particular, 
CO"*" is mainly produced via the reaction C+-I-OH— >^CO''"-|-H. The production of OH is very 
dependent on the temperature as O+H2 may come into play at the H/H2 transition region 
when H2 is abundant and the temperature is still a few hundred K. The corresponding 
endothermicity is about 3000 K. An increase in the gas kinetic temperature in the HI/H2 
interface could have a dramatic effect in the CO"*" production. 

There are several observational evidences that suggest that chemical models fail to 
predict the gas kinetic temperature in the HI/H2 region. One of the best studied PDRs 
is the Orion Bar, which is the paradigm of a Galactic PDR associated to an HIl region. 
Observations of the H2 rotational lines in the Orion Bar by Parmar et al. (1991) revealed 
the existence of unexpectedly large amounts of warm gas (T~400-700 K). They proposed 
that dumpiness could help to reconcile the observations with chemical models. Recently, 
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AUers et al. (2005) proposed that the dust FUV attenuation cross sections should be reduced 
by a factor of 3 in order to explain the separation between the ionization front and the H2 
emission peak in the Orion Bar. In order to explain the intensities of the H2 ro- vibrational 
lines they need to re-adjust the photoelectric heating rate. A change in the size distribution 
of the grains and/or the photoelectric heating rate would produce large variations in the 
thermal balance of the PDR. 

This paper has been partially funded by the Spanish MCyT under projects DGES/AYA2000- 
927, ESP2001-4519-PE, ESP2002-01693, AYA2002-01241 and AYA2003-06473. 
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Fig. 1. — Observed spectra towards the positions East (E), Center and West (W) in M82. 
The 30m beam at 236 GHz around the three observed positions has been drawn in the 
interferometric H^^CO"'" image by Garcia-Burillo et ah (2002). 
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Fig. 2. — Model predictions for the fractional abundances (a.)and cumulative column den- 
sities of CO"*" and HCO"^ (b.) derived using the Le Bourlot et al.'s code for the physical 
conditions in the M 82 nucleus (Go=10^ in units of the Habing field, n=10^ cm~^ and C=4 
10-^^ s"^). Note that the N(CO+)/N(HCO+) ratios measured in the M82 nucleus (horizon- 
tal dashed lines in c. panel) are well explained if the emission arises in PDRs with total 
visual extinction between 4.5 and 6.5 mag, in perfect agreement with our results in Paper I 
based on the CN/HCN ratio. 
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Table 1: Observational parameters and Gaussian fits results 





Rest Frcq." 


/ TMsdv 


Vlsr 


A?; 




t'tnt 


Molecule 


(MHz) 


(K km s"^) 


(km s~^) 


(km s~^) 


(mK) 


(min) 


E. (+14", +5") 


C0+ J=5/2-^3/2 


236062.55 


0.36(0.08) 


311(8) 


74(19) 


4.6 


695 


HCO+ J=3^2 


267557.00 


31.33(0.45) 


301(1) 


96(2) 


304.7 


127 


HOC+ J=3^2 


268451.00 


0.65(0.13) 


307(5) 


51(12) 


12.0 


125 


fji3co+» J=1^0 


86754.33 


2.05(0.15) 


320(2) 


66(6) 


29.0 




Center (0",0") 



C0+ J=5/2- 


-^3/2 


236062.55 


0.88(0.11) 


181(9) 


137(19) 


6.0 


340 


C0+ J=3/2- 


-.1/2 


235789.64 


0.34(0.09) 




137^ 


2.4 


340 


HCO+ J=3- 


-^2 


267557.00 


34.66(1.52) 


147(4) 


149(7) 


218.1 


12 


HOC+ J=3- 


-^2 


268451.00 




<0.34^ K km 






190 


gi3(30+« J= 


= 1-^0 


86754.33 


1.21(0.13) 


171(8) 


140(18) 


8.1 




CH3OH Jk= 


52,3~>'42,2'^ 


241904.63 


0.92( 0.27) 


222(13) 


85(27) 


10.0 


70 



W. (-14" ,-5") 


C0+ J=5/2^3/2 
fji3co+» J=l^0 


236062.55 0.38(0.10) 
86754.33 1.94(0.15) 


111(8) 
103(3) 


81(21) 4.4 390 
73(7) 24.8 



" Spectra obtained from the interferometric H^'^CO^ image reported by (Garcfa-Burillo et 
al. 2001) by convolving to a Gaussian resolution of HPBW=11". 
''Gaussian fit obtained by fixing the velocity and linewidth. 
'^3x(j limit with Av=50 km s~^ 
The fine Jfc=5-2,4^4_2,3 at 241904.15 MHz is overlapped with this. 



Table 2: Column densities and relative fractional abundances 



Molecule 




M82i 




Orion Bar^ 


NGC 7023^ 


Mon R2'* 




E. 


(0,0) 


W. 


IF 


PDR peak 


IF 


N(CN) 


6.3 xlO^""' 


8.8x10^5 


l.lxlO^^* 


2.4x10^-* 


2.4x10^-* 




N(CO+) 


l.SxlQi^ 


3.7xl0i3 


1.6x10" 


8.0x10" 


1.6x1012 


4.4xl0i2 


N(HOC+) 


2.5x10" 


4.3x10" 


2.5x10" 


4.0x10" 


1.3x10" 


1.4xl0i2 


CN/HCN 


6 


8 


7 


3 


4-8 




HCO+/HOC+ 


44 


36 


30 


<166 


50-120 


460 


CO+/HCO+ 


0.04 


0.1 


0.04 


~0.01 


0.01-0.11 


0.005 



1 Assuming a size of 6" and data from Paper I and this paper. 

2 Assuming that the bar fills half of the beam and data from Fuente et al. (1993,1996b). 
^ Assuming a filament of 6" and data from Fuente et al. (1993,2003). 

* Assuming the size of the radio continuum emission at 6 cm, 10", and data from Rizzo et al.(2003). 



